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Abstract 

We discuss the phenomenology of damping signatures in the neutrino oscillation probabili- 
ties, where either the oscillating terms or the probabilities can be damped. This approach 
is a possibility for tests of non-oscillation effects in future neutrino oscillation experiments, 
where we mainly focus on reactor and long-baseline experiments. We extensively motivate 
different damping signatures due to small corrections by neutrino decoherence, neutrino 
decay, oscillations into sterile neutrinos, or other mechanisms, and classify these signatures 
according to their energy (spectral) dependencies. We demonstrate, at the example of short 
baseline reactor experiments, that damping can severely alter the interpretation of results, 
e.g., it could fake a value of sin 2 (29 13) smaller than the one provided by Nature. In ad- 
dition, we demonstrate how a neutrino factory could constrain different damping models 
with emphasis on how these different models could be distinguished, i.e., how easily the 
actual non-oscillation effects could be identified. We find that the damping models cluster 
in different categories, which can be much better distinguished from each other than models 
within the same cluster. 
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1 Introduction 



Neutrino oscillations are by far the most plausible description of transitions among different 
neutrino flavor eigenstates [1-7]. However, there have historically been other attempts 
in the literature to describe these transitions with other mechanisms as well as neutrino 
oscillations combined with such other mechanisms. These scenarios include neutrino wave 
packet decoherence [8-12], neutrino decay [13-20], oscillations into sterile neutrinos [21,22], 
neutrino absorption (see, e.g., Ref. [23]), and neutrino quantum decoherence [24-32]. A 
combined scenario is, for example, the combination of neutrino oscillations and neutrino 
decay (see, e.g., Refs. [19,20]). Although these other mechanisms, leading to "non-standard 
effects", are not such successful descriptions for flavor transitions as neutrino oscillations 
are (in fact, they are strongly disfavored [6,7]), they could still give rise to small corrections 
to the neutrino oscillations. These non-standard effects need to be described in a framework 
together with neutrino oscillations and can be constrained by current and future experiments 
(see, e.g., Ref. [33] for a recent review). Thus, we will assume that the leading order effect in 
neutrino flavor transitions is due to neutrino oscillations, whereas the next-to-leading order 
effects are described by different "damping mechanisms" of the neutrino oscillations. 

Since any non-standard effect may point towards new interesting physics beyond the stan- 
dard model, the test of small corrections due to these effects should be one of the main 
objectives in future high-precision neutrino oscillation physics. The assumption of stan- 
dard three-flavor neutrino oscillations will inevitably lead to an erroneous derivation of the 
elements of the mixing matrix U or the mass squared differences. We therefore define 
"non-oscillation effects" as any modification of the three-flavor neutrino oscillation proba- 
bilities in vacuum as well as in matter. For example, the LSND anomaly [34] could be an 
indication of non-oscillation effects according to this definition. Since future reactor and 
long-baseline neutrino oscillation experiments are expected to have high-precisions to the 
subleading neutrino oscillation parameters sin 2 (26> 13 ) and <5cp, we mainly discuss the impact 
of non-oscillation effects or possible constraints on the non-oscillation effects in the context 
of these experiments. 

In principle, one could think of many different approaches to test non-oscillation effects with 
future long-baseline experiments: 

Neutral-currents can be used to test the conservation of probability, i.e., P ae + P a[l + 
P aT = 1 (see, e.g., Ref. [35]). However, at long-baseline experiments, uncertainties 
in the neutral-current cross-sections and the charged-current contamination lead to a 
precision of only about 10 % — 15 % [35]. In addition, even if some non-oscillation 
effects are found, there will be no information on the nature of the effects, whereas 
effects conserving the overall probability cannot be detected at all. 

The detection of v T can complement the information on P ae and P afl to test the conser- 
vation of probability (see, e.g., Ref. [36]). Since v T detection is much more sophisti- 
cated and less efficient than the detection of v e and due to the higher r production 
threshold, this is also a non-trivial test. If there are non-oscillation effects, then the in- 
formation will be better than in the preceding case, since one will know which neutrino 
oscillation probabilities are affected. 
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Unitarity triangles for the lepton sector can be constructed [37, 38] . However, since 
there is no simple relationship between the quantities of the unitarity triangles and 
the neutrino oscillation observables, this approach may not be the most feasible for 
the lepton sector. 

Tests of distinctive signatures, i.e., spectral (energy) dependent effects, could directly 
identify certain classes of non-oscillation effects [33,39-41]. The advantage of such 
tests is that the effect could be directly identified if it produces a unique signature in 
the energy spectrum. In addition, this test does not depend upon normalization errors 
of the event rates, which are likely to constrain the first two measurements. However, 
there might be strong correlations with the neutrino oscillation parameters. 

In addition, in the future, it may be possible to resolve the line width and shape of the 7 Be 
solar neutrino line [42,43] and extract the temperature distribution as well as the modulation 
of this line, which could be caused by next-to-leading order effects. Thus, performing very 
high-energy resolution measurements of the 7 Be line may be an idea how to determine 
these next-to-leading order effects. Such possible precision neutrino experiments include, 
for example, a bromine cryogenic thermal detector proposed in Refs. [44,45]. 

In this study, we will focus on the tests of distinctive signatures in which we introduce 
"damping signatures" as an abstract concept for a class of possible effects entering at prob- 
ability level. 1 In general, small Hamiltonian effects, see, e.g., Ref. [28], may be as important 
as the kind of damping effects that we will describe in this study. Such Hamiltonian effects 
could lead to direct changes in the effective neutrino oscillation parameters. Nevertheless, 
those effects cannot be treated in the framework presented here. We will use the observation 
that mechanisms, such as decoherence or decay, lead to exponential damping in the neutrino 
oscillation probabilities. However, the effect might be stronger for low or high energies, i.e., 
the spectral (energy) dependence for the damping might be different. A common feature 
for many of the discussed models is that they will lead to less neutrinos (of all active fla- 
vors) being detected than what is expected from the three-flavor neutrino oscillations. For 
all other models, only the oscillating terms of the neutrino oscillation probabilities will be 
damped, while the total number of active neutrinos remains constant. Note that the damp- 
ing signature approach does not cover all possible models, but many models can, at least 
in the limit of small corrections, lead to some exponential damping effect. 

Our study is organized as follows. In Sec. El we will present and classify different forms 
of the damping signatures. For the reader who is not interested in different models for 
damping signatures, at least Sec. 12.11 and the examples in Table [T] should be read to be 
able to follow the rest of the study. Next, in Sec. El we will give and discuss the damped 
neutrino oscillation probabilities arising from the effects described by their signatures. For 
the reader, who is most interested in possible experimental implications, Sec. 13.11 should 
summarize the most relevant features, whereas the rest of this section deals with the more 
technical three-flavor cases. Then, in Sees. H] and EJ we will discuss the physics of these 
damping signatures and give two different applications in the framework of a complete 
experiment simulation. Especially, in Sec. 01 we demonstrate how such damping signatures 

1 Although it will be possible to describe some of our effects on Hamiltonian level, the Hamiltonian will 
not be Hermitian anymore. 
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can modify the interpretation of physical results for future reactor experiments, whereas in 
Sec. we discuss how a neutrino factory could constrain different damping signatures and 
how these different signatures could be distinguished. Finally, in Sec. El we will summarize 
our work and present our conclusions. 



2 Phenomenology of damping signatures 



In this section, we motivate, in a phenomenological manner, the form of the damping sig- 
natures used for the rest of this study. 



2.1 General description of damped neutrino oscillations in vacuum 



We start with three-flavor neutrino oscillations in vacuum, which can be described by the 
(undamped) vacuum oscillation probabilities 



IMU diag (l,exp (-i^) ,exp (-i^)) ^> 

3 

= J2 U ^ U h U * ai UpieM-^ij)- (1) 

Here U is the leptonic mixing matrix in vacuum, Am| = m| — m 2 - the mass squared 
difference, and $jj = AmfjL/(2E) the oscillation phase. By defining 

n R * * » Am%L mf - m 2 

Jg ^UaiUyrjJf* and ^ = -^ = -L-—L L = -± t 

the oscillation probabilities may be written as 

3 

Pap = ^Re(J^)-4 Re(J^)sin 2 (A,,)-2 £ Im( J$>) sin(2A i ,) 

ij'=l l<«<j<3 l<i<j<3 



$>f + 2 £ |J<f|cos(2A^ + argj; 

i=l 1<«<J<3 



ij ) 



where, in the first line of the equation, the first two terms are CP-conserving and the third 
term is the source of any CP violation, this corresponds to arg being the source of any 
CP violation in the second line. As will be discussed, there may be reasons to assume that 
Eq. (J2J) does not give the correct neutrino oscillation probabilities. Effects that might spoil 
this approach of calculating neutrino oscillations probabilities include loss of wave packet 
coherence and neutrino decay. The effective result of such processes is to introduce damping 
factors to the oscillating terms of the neutrino oscillation probabilities. We define a general 
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damping effect to be an effect that alters the neutrino oscillation probabilities to the form 



P af 3 = Yl Ua i U h U ai U Pi exp(-i^) 



3 



^Re(4f)A,--4 Re(J^)A, sin 2 (A^-2 £ Im(jf )D ij sin(2Ay) 

ij'=l l<i<i<3 1<«<J<3 



3 



^jfAi + 2 £ |J^|AiCos(2A, i + argJ^) ) (3) 

i=l l<«<j<3 



where the damping factors 



Aj = exp [ -a^ -jL ) (4) 



have been introduced and we have assumed that Aj = Dji- Obviously, as Aj — > 1, we 
regain the undamped oscillation probabilities given in Eq. (J2J). In Eq. (j3J), > is a 
non-negative damping coefficient matrix, and (3, 7, and £ are numbers that describe the 
"signature", i.e., the L and E (7) dependencies as well as the dependence on the mass 
squared differences. In addition, the parameter £ implies two interesting cases: 

£ > 0: In this case, only the oscillating terms will be damped, since Am^ = by definition. 

£ = 0: The whole oscillation probability can be damped (depending on Oiij), since also the 
terms which are independent of the oscillation phases are affected. 

Therefore, we expect two completely different results for these two cases. In general, Eq. (j3J) 
introduces twelve new parameters, which can be used to model many non-standard contri- 
butions that enter on the oscillation probability (not Hamiltonian) level. We will give some 
examples of such contributions below. Although we expect these contributions to be small, 
it is rather impractical to deal with that many new parameters, which means that some 
simplifications need to be made. First of all, note that the parameter (3 is not measurable 
if only one baseline is considered and can therefore be absorbed in atij. For two baselines, 
it can, in principle, be resolved if all the other parameters are known. Second, for a specific 
model, there may be relations among different ay's that actually imply much fewer inde- 
pendent parameters. For a very simple model, the number of parameters can even reduce 
to one. Since we are mainly interested in the spectral signatures, i.e., 7, we will often use 
otij = a to estimate the magnitude of different effects. Third, it will turn out that the 
parameter £ is strongly dependent on the model, since, as discussed above, it describes two 
completely different classes of models. Hence, we will finally end up with one free parameter 
a and several fixed model dependent parameters (3, 7, and £. 

2.2 A model for damped neutrino oscillations in matter 

In some cases, we will use neutrino propagation in matter, since, for instance, neutrino 
factories operate at very long-baselines for which matter effects become important. We 
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use an approach similar to Eq. (J3J), which should describe the damping signatures as minor 
perturbations to neutrino oscillations in (constant) matter as long as they are small enough: 

3 

P aP = £ Re(JJf) A, - 4 Re(J^)A,sin 2 (A J ,)-2 £ Im( jf ) D y sin(2A y ), 

ij=l l<j<?'<3 l<«<i<3 

(5) 

where the tildes denote the effective parameters for neutrinos propagating in matter (for 
instance, J°f = UajU^U^Upi, where U is the effective leptonic mixing matrix in mat- 
ter, i.e., the matrix re-diagonalizing the Hamiltonian with the matter potential included). 
In general, the damping effects may not enter directly as multiplicative factors in the in- 
terference terms among different matter eigenstates. 2 However, in this study, we assume 
small damping effects that should act as perturbations which, to leading order, give rise to 
neutrino oscillation probabilities in matter of the same form as the ones in vacuum. 

Thus, we use the propagation in constant matter and apply the damping signatures to the 
mass eigenstates in matter. This means that we discuss signatures which depend on the 
mass eigenstates in matter. They may come from wave packet decoherence, neutrino decay, 
neutrino oscillations into sterile neutrinos, neutrino absorption, quantum decoherence, or 
other mechanisms. Strictly speaking, this model does not describe many of these mecha- 
nisms exactly, since a complete re-diagonalization of the Hamiltonian might be necessary 
(such as for Majoron decay in matter; see, e.g., Refs. [46,47]). However, we treat only small 
effects in matter acting as a perturbation to the neutrino oscillation mechanism and do not 
consider transitions from active into active neutrinos, which would require a more compli- 
cated treatment (such as decay into other active neutrino states). Therefore, this model 
should be sufficient as a first approximation, since we will later on use either short baselines 
or mainly discuss effects in the P w channel, which are not affected by matter effects to first 
order in the ratio of the mass squared differences Am^/Arojj and the mixing parameter 
s 13 = sm(9 13 ) [48]. 

2.3 Examples of different damping signatures 

The general damping signature in Eq. (j3J) seems to be very abstract. Therefore, let us 
now give some motivations for such damping signatures by different mechanisms, which are 
summarized in Tabled 

Intrinsic wave packet decoherence 

Intrinsic wave packet decoherence is an effect that appears even in standard neutrino os- 
cillation treatments [8-12]. It naturally emerges from any quantum mechanical model that 
does not assume neutrino mass eigenstates propagating as plane waves or from any quantum 
field theoretical treatment. In principle, intrinsic decoherence may not be distinguishable 
from a macroscopic energy averaging (see, e.g., discussions in Refs. [49-51]). Therefore, it 

2 For instance, some effect on Hamiltonian level, such as neutrino absorption, would require a full re- 
diagonalization of the effective Hamiltonian with the absorption terms included, see the section "Neutrino 
absorption" below. 
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Damping type Signature D 



VI 



Unit for a 


(3 


7 




MeV 2 or GeV 2 


2 


4 


2 


GeV ■ km" 1 


1 


1 





eV 4 


2 


2 





GeV" 1 • km" 1 


1 


-1 





GeV" 2 • km" 1 


1 


-2 





eV" 2 


1 or 2 


2 


2 



Wave packet / 2 (Am|) J i ! 
decoherence 6X ^ V ^ 8£4 

Decay exp (— a|u 

Oscillations to v s exp y—C j^i 

Absorption exp (—aLE) 

Quantum , 2 , 

, , T exp ( — aLh/ 
decoherence 1 

Quantum / (Am^ 

decoherence II eX ^ \ K E ' 2 



Table 1 : Different examples for damping signatures considered in this study. The parameter 7 represents 
the spectral (energy) dependence of the signature. The parameter a has in some places been re-defined for 
convenience (see main text) unless it corresponds exactly to our definition of a. The quantum decoherence 
models I and II are two examples of signatures motivated by quantum decoherence (see Table |2J ■ The 
quantum decoherence model II absorbs (3 in the definition of k = aL 13 in order to describe two of the 
models from Tabled Note that another commonly used quantum decoherence signature is the same as the 
decay signature. 



is natural to expect that the test of this signature could be limited by the knowledge on the 
energy resolution of the detector. 

We adopt the treatment in Ref. [8], which uses averaging over Gaussian wave packets. In 
this approach, the loss of coherence can only be described at probability level. It leads to 
factors exp [-(L/L^ h ) 2 ] in Eq. (gj), where L-° h = 4y/2a x E 2 /\ Am 2 .| and a x is the spatial 
wave packet width. In this case, the damping descriptions in vacuum and matter using 
Eqs. (JSJ), 0? and © are accurate. For the damping signature, we obtain 



D 



ij 



exp 













= exp 







\[2a E Am%L 



E 



AE 



exp 



-4 



(Am; 



2 \2 7-2 
ij) L 



8£ 4 



(6) 



in vacuum and the analogous signature D in matter. Here we have introduced a wave packet 
spread in energy <je = l/(2a x ), since we later will derive an upper bound for this quantity 
and directly compare it to the energy resolution of a detector. The typical units of ge will 
be MeV or GeV. By comparing Eqs. (HJ) and ©, we can identify a^- = cr%/8, (3 = 2, 7 = 4, 
and £ = 2. Note that, in this case, the cWs do not depend on the indices i and j. 



In order to better understand Eq. 
first oscillation maximum: 



we note that Am 2 L/(4£') is of order unity for the 



D 



ij 



exp 



f V2a E AmfjL 



E 



AE 



exp 



V2E 3 



exp 



V2a x E 



0(1; 



value at oscillation maximum 



(7) 

From Eq. (|7J), we find three major implications: First, it means that no effect will be observed 
if oe "C E, because the oscillation phase is usually of order unity (or less). Second, since the 



6 



decoherence damping factor always comes together with an oscillation phase factor with the 
same Am^ [cf., Eq. ©], it will equally damp the solar and atmospheric oscillating terms in 
one probability formula. This means for the atmospheric oscillation experiments that if the 
solar contribution cannot be neglected, its damping factor can also not be neglected. Third, 
one expects the largest suppression for low energies independent of the type of oscillation 
experiment (solar or atmospheric), since in either case the experiment will be operated 
close to the oscillation maximum. Eventually, it is important to keep in mind that this 
decoherence signature is not an intrinsic property of the neutrinos, but an effect related to 
the production and detection processes. Therefore, the parameter oe could be different for 
different classes of experiments. 

Invisible neutrino decay 

Another example of a damping signature is neutrino decay (see, e.g., Refs. [13-18]). In 
particular, invisible decay, i.e., decay into particles invisible for the detector, leads to a 
loss of three-flavor unitarity. In this case, the neutrino evolution is given by an effective 
Hamiltonian 



where Y = diag(ai, a 2 , a 3 )/2 in the neutrino mass eigenstate basis, = Ti/ji, Tj is the 
inverse life-time of a neutrino of mass eigenstate i in its own rest frame, and 7« = E/rrii 
is the time dilation factor. We note that H and T are both diagonal in the neutrino 
mass eigenstate basis. The neutrino oscillation probabilities may now be calculated as 
usual with the exception that, in addition to the phase factor exp[— \m\Lj {2E)\, a factor of 
exp[— TirriiL/ (2E)] is obtained when evolving the neutrino mass eigenstate z/j. The resulting 
neutrino oscillation probabilities are of the form of Eq. (jHJ) with 



where a.i = Tirrii, in accordance with Refs. [19,20]. Thus, for neutrino decay, the character- 
istic signature is a^- = (a« + atj)/2, ft — 7 = 1, and £ = 0. 

An example of the above decay is Majoron decay into lighter sterile neutrinos. In this case, it 
is plausible to assume a quasi-degenerate neutrino mass scheme for the active neutrinos with 
approximately equal decay rates for all mass eigenstates, since the decay products all have 
to be considerably lighter than the active neutrinos to obtain fast decay rates due to phase 
space. The decay rates of the cVs will then be approximately equal (a; = a for all i) and 
will typically be given in units of GeV/km. Note that the decay rate is an intrinsic neutrino 
property, not an experiment-dependent quantity such as the wave packet decoherence. We 
identify by the comparison of Eq. (0) with Eq. (JU) that a is the same quantity 3 , ft = 7 = 1, 
and £ = 0. In matter, we use the analogous signature, i.e., we let the mass eigenstates in 
matter decay. In general, this is only a first approximation, since, for example for Majoron 
decay in matter, a re-diagonalization of the complete Hamiltonian may be necessary; see, 
e.g., Refs. [46,47]. However, as we have assumed equal decay rates for all eigenstates, it 

3 In general, we do not change the symbol for a if its is exactly the same as the one in Eq. Q. However, if 
there are additional factors absorbed in a, then we re-define the name (such as for wave packet decoherence). 



H C B — H — iV, 





(9) 
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should describe the problem exactly, since the mass eigenstates in matter will also decay 
with equal rates. In different decay models, the a^'s may not be identical anymore. For 
example, for a hierarchical mass scheme with a normal hierarchy, the mass eigenstate m 3 
decays much faster than the other two. In this case, the observed effects in atmospheric 
oscillations would qualitatively be similar, but about a factor of two smaller (since mainly rri2 
and 7713 participate in the oscillation and only one of them decays). However, in matter such 
a model is much more difficult to treat, since it is not easy to identify the mass eigenstate 
in matter after the diagonalization of the Hamiltonian. This problem does not occur with 
equal decay rates. 



Oscillations into sterile neutrinos 

A natural description for the LSND result [34] is a light sterile neutrino (i.e., not a weakly 
interacting neutrino) that is mixing with the active neutrinos. This description is now dis- 
favored for the LSND experiment [21,22], but small admixtures of light sterile neutrinos 
cannot be entirely excluded. In particular for slow enough oscillations into sterile neu- 
trinos, the oscillation signature sin 2 A^ with Ay = Am?jL/ (AE) translates into damping 
signatures: 



1 — e sin^ 



Am|L 
AE 



1 - e 



Am|L 
AE 



cxp 



Am|L 
AE 



(10) 



where e represents the magnitude of the mixing. Thus, the damping coefficient a will (in 
this case) be determined by the sizes of the mixing and the mass squared differences Am|. 
We use as a model in vacuum (and the same form in matter) 



Dij = exp [ -a Vi0 



(2Ef 



exp 



(2Ef 



(11) 



where e contains the information on mixing and Am 2 and will be given in units of eV 4 (the 
mixing factor is dimensionless) . Thus, we identify by comparison of Eq. (|11|) with Eq. P|) 
that ctij — e/4, (3 = 7 = 2, and £ = 0. Note that we only discuss effects independent of i 
and j, which simplifies the problem, but restricts the number of applications tremendously. 
In addition, although the coefficient e is not experiment dependent (since it is an intrinsic 
neutrino property here), it may (partly because of the independence on i and j) depend on 
the oscillation channel and mass scheme. As an example, let us consider P w and a mass 
scheme with Am 2 x <C Am 2 ^ < Am 2 ^ i.e., Am^ is the largest mass squared difference. In 
this case, one can show that to first approximation e ~ f/ 2 4 U^ 3 (Aml 3 ) 2 (for CP conserva- 
tion). Thus, e is suppressed by the flavor content of in and the extra mass squared 
difference, since all the other mass squared differences with the sterile state are absorbed 
into the atmospheric oscillation terms. In general, it should be noted that sterile neutrinos 
are not affected in the same way as active neutrinos when propagating through matter (i.e., 
there is a phase difference due to the neutral-current interactions between matter and the 
active neutrino flavors). However, the exponential damping signature for oscillations into 
sterile neutrinos presented here is only valid for short baselines, where matter effects have 
not yet developed. 



8 



Neutrino absorption 



When neutrinos propagate through matter, there is a small chance of absorption. Neutrino 
absorption can be described in a fashion similar to neutrino decay. In this case, we assume 
that an effective Hamiltonian is given by 

H cS = H- ir, (12) 

where H is the usual neutrino Hamiltonian in matter, V is given by 

r = pdiag(a e ,cr fl ,a T )/2 (13) 

in the flavor eigenstate basis, p is the matter density, and a a is the absorption cross-section 
for a neutrino of flavor a. If we assume the cross-sections to be relatively small, then the 
eigenstates of H e g will not differ significantly from the orthogonal eigenstates of H . Thus, 
the first order corrections to the eigenvalues of the effective Hamiltonian will be 

SEP = -iT u = -i| Yl \ U <*i\ 2 <?a = -fyi, (14) 

a 

where <j{ is an effective cross-section for a neutrino of mass eigenstate i. The neutrino 
oscillation probability is now given by an expression of the form of Eq. (J3J) with 

•m±^V) , (is, 

where we have assumed a constant matter density p. The signature of this scenario is given 
by (3 = 1 and 7 is equal to minus the power of the energy dependence of the cross-sections. 
It should be observed that, since the cross-sections increase with energy, 7 will be a negative 
number. 

If all neutrino flavor cross-sections were equal (or approximately equal), then the effective 
matter eigenstate cross-sections would also be equal. 4 For the neutrino energies relevant to a 
neutrino factory, the neutrino-nucleon cross-sections are approximately linear in energy [55]. 
Thus, in this energy range, the damping signature is given by a = pa(E )/E , (3 = 1, 
7 = —1, and £ = 0, where a(E ) is the cross-section at energy E . At higher energies, the 
cross-sections increase at a slower rate and if damping effects are studied at these energies, 
then the effective damping parameter 7 lies in the interval — 1 < 7 < 0. 

It should be noted that the standard neutrino absorption effects (by weak interactions) are 
very small for energies typical for neutrino oscillation experiments. However, there could 
be non-standard absorption effects and the cross-sections of these effects should behave in 
a manner similar to the standard absorption. 

4 Because of the higher r production threshold, the v e and cross-sections are in fact considerably 
larger than the v T cross-section [52-54]. However, for these low energies the standard absorption effects are 
anyway small. 



D 



exp 



&i + o-j 



pL 



exp 
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Reference Signature Dij Unit for a (3 7 £ 

Lisi et al. [24] and Morgan et al. [32] exp (-aL) km -1 1 0~ 

Lisi et al. [24] and Morgan et al. [32] exp (-afi) GeV • knT 1 1 1 

Lisi et al. [24] and Morgan et al. [32] exp (—aLE 2 ) GeV -2 • km" 1 1 -2 

Adler [26] exp [~a {Am ^ L ^ GeV" 1 1 2 2 

/ (Am 2 ) 2 L 2 \ 

Ohlsson [27] exp ( — ex — jjj^ j dimensionless 2 2 2 

Table 2: Different signatures that might arise from quantum decoherence and the references in which 
they are motivated. 



Quantum decoherence 

It has been argued that quantum decoherence could be an alternative description of neutrino 
flavor transitions. Fits to data by different collaborations (e.g., Super-Kamiokande [7] and 
KamLAND [6]) have been performed and these clearly disfavor a decoherence explanation 
for neutrino flavor transitions. However, quantum decoherence may still be a marginal effect 
in addition to neutrino oscillations and could give rise to damping factors of the type given 
in Eq. (@J). 

Quantum decoherence arises when a neutrino system is coupled to an environment (or a 
reservoir or a bath), which could consist of, for example, a space-time "foam" [24] leading 
to new physics beyond the standard model. Thus, quantum decoherence may be a feature 
of quantum gravity. In order to find the formulas describing quantum decoherence, it is 
necessary to use the Liouville equation with decoherence effects of the Lindblad form [56] . 

Throughout the literature [24-32,57,58], the effects of loss of quantum coherence in neutrino 
oscillations have been studied. Although the signatures derived by different authors seem 
to vary, the decoherence effects are of the same form as Eq. (J1J). However, there might be 
additional effects on the oscillation phases. In Table El we give a brief summary of some of 
the signatures that are present in the literature, these examples could be used to motivate 
the numerical testing of such signatures. 



Other signatures 

In principle, what we have presented above is just a collection of interesting signatures 
that could be responsible for damping of neutrino oscillations. However, there are also 
other possibilities, which we have decided not to investigate further in this study. These 
signatures include, for example, heavy isosinglet neutrinos [59, 60] and neutrino oscillations 
in different extra dimension scenarios [32,61-65]. 



Combined signatures 

In most cases, if there is a damping effect, then it would be natural (and easy) to assume 
that one type of effect is giving a clearly dominating contribution. However, if an experiment 
is carried out with some specific setup, then contributions from different scenarios might be 
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of the same order. In such a case, the form of Eq. (j3J) is spoiled. For example, in the case 
of neutrino decay combined with neutrino absorption, the matrices T are just added which 
results in the damping signatures 

Dij = exp 

In general, just multiplying the damping factors (which is the result of the above treatment) 
might not give the correct damping and different combined cases might behave in other ways. 
However, since there are different energy dependencies in the different damping signatures, 
there will only be a limited energy range where a combined treatment is necessary. In this 
study, we do not consider combined signatures. 

3 Damped neutrino oscillation probabilities 

In this section, we investigate the effects of damping on specific neutrino oscillation prob- 
abilities interesting for future reactor and long-baseline experiments, where we restrict the 
analytical discussion to the vacuum case. 

3.1 The damped two-flavor neutrino scenario 

In a simple two-flavor scenario, the damped neutrino oscillation probabilities take partic- 
ularly simple forms (just as in the non-damped case). From the two-flavor equivalent of 
Eq. (jSJ), we obtain 

P aa = D u c A + D 22 s 4 + ^D 21 sin 2 (20) cos(2 A), (17) 

Ppp = Dns 4 + D 22 c 4 + ^D 21 sin 2 (26) cos(2 A) (18) 
for the neutrino survival probabilities and 

P a /3 = Pp* = \ sin 2 (26) [D n + D 22 - 2D 21 cos(2A)] (19) 

for the neutrino transition probability, where v a is the linear combination v a = cv\ + su 2 , 
Up is the linear combination that is orthogonal to u a , A = A 2 i, s = sin(#), c = cos(#), and 
6 is the mixing angle between the two neutrino flavors. 

Let us first discuss the case £ > or all an = 0, which means that all Da are equal to 
unity. We refer to this case as "decoherence-like" (probability conserving) damping. The 
two-flavor formulas then become 

Pa/3 = S af 5 + \(l ~ 2M sin 2 (2#)[l - £>cos(2A)], (20) 

where D = D 2 i. Below, we will show that expressions reminding of these two-flavor formulas 
will be quite common in the three-flavor counterparts. In the limit D — » (maximal 
damping), the oscillations are averaged out, i.e., 

Pa/3 -> Ml - sin 2 (26)} + i sin 2 (26), 



a 



decay 



+ a^E I L 



E 



(16) 
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where the factor 1/2 is typical for an averaged sin 2 (x) term. It is also of interest to note, 
from the form of Eq. ([20jh that the neutrino transition probabilities can either be smaller 
or larger than the undamped probabilities depending on the sign of cos(2A). For instance, 
the neutrino survival probability 



1 - -sin 2 (2#)[l - Dcos(2A)]. 



(21) 



is smaller than the corresponding undamped probability if cos(2A) is positive and vice 
versa. Close to the oscillation maximum A ~ 7r/2, the factor cos(2A) will be negative, i.e., 
the damped neutrino survival probability will be larger than the undamped probability, 
since the oscillations will be partially averaged out. This behavior chang function 
of the neutrino energy at points where cos(2A) changes sign, i.e., at 2A = (n + l)7r/2, 
n = 0, 1, . . .. As a rule of thumb, the damping will lead to larger probabilities close to the 
oscillation maximum E max = Am 2 L / (2n) and to smaller probabilities for E < 2E max /3 and 
E > 2E ma , x . This result will be valid for any survival probability discussed in this study. 

From the form of Eq. ()20|) . it is apparent that if only a small range of A's is studied, then 
a damping factor may mimic an oscillation signal. The worst such case would be if the 
damping signature had 7 = 2. This would mean that if one makes a series expansion of 
cos(2A) and the exponential of the damping factor, then the energy dependence will be the 
same to lowest order in the expansion parameters, i.e., we will have 



Dcos(2A) 



,A 2,^ 

1 -a\Am 2 f — + 



1 



Am 2 L 
AE 



+ 



(22) 



This effect is also present in a general case with any number of neutrino flavors. 

Another interesting case is when = a, + <x,- and £ = 0, which is expected for the 
neutrino decay and neutrino absorption scenarios. This assumption results in the fact that 
the damping factor Dij can be written product 



D 



1.1 



AiAj, 



(23) 



where A4 = exp(— aiL 13 / E 1 ) is only dependent on the ith. mass eigenstate. Then, the 
neutrino oscillation probabilities are given by 



p 

1 aa 

Pp/3 

Pa(3 



A 2 [{c 2 + Ks 2 ) 2 - /tsin 2 (2#) sin 2 (A)] 
A 2 [(kc 2 + s 2 ) 2 - Ksin 2 (2#) sin 2 (A)] 

- A 2 sin 2 (26) [1 + k 2 - 2k cos(2 A)] , 



(24) 
(25) 

(26) 



where A = A\ and k = A^j A\. It is important to note that, for example, the total 
probability P aa + P a p is not conserved in this case, in fact, we obtain 



E*aa Paf3 



A 



c 4 + K V + ^sin 2 (20)(l + K ) 2 



< 1. 



(27) 



where the equality holds if and only if A — k — 1 (because of the form of the A^s, A < 1, 
kA < 1, and that all terms in Eq. ()27|) are positive, the terms will attain their maximum 
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value when A = kA = 1, in which case the entire expression simplifies to one). Thus, we 
will introduce the term "decay-like" for effects giving rise to damping terms of the form 
given in Eq. ([23~j). 

In the case of a decay-like signature, there are two special cases which are of particular 
interest. First, if both mass eigenstates are affected in the same way, i.e., k — 1, then 
the resulting neutrino transition probabilities will reduce to the undamped standard neu- 
trino oscillation probabilities suppressed by a factor of A 2 . This means that all damped 
probabilities will be smaller than their undamped counterparts. Second, if only one of the 
mass eigenstates is affected, i.e., A = 1, then the difference in the v a survival probability 
compared to the undamped case will be given by 

AP QQ = P Q t mped - P^ damped = (k - l)s 2 [(l + k)s 2 + 2c 2 cos(2A)]. (28) 

Thus, this survival probability will actually increase if 

-2 cos(2A) > (1 + k) tan 2 (#). (29) 

Note that for the first part of the neutrino propagation (for L < 7cE/Am 2 ), the term cos(2A) 
is positive, and thus, the inequality of Eq. (j2Hj) cannot be satisfied in this region, since the 
right-hand side is always positive. From the comparison with the discussion after Eq. (j21j) . 
this condition is equivalent to E > 2P max . For example, for a neutrino factory, which can 
be operated far away from the oscillation maximum, this implies that the relevant part 
of the spectrum will be suppressed by this form of damping. For the neutrino oscillation 
probability difference AP Q/ g, we obtain 

AP a/3 = isin 2 (2#)(K-l)[l + fi;-2cos(2A)], (30) 

that is, the damped P a/ 3 is larger than the undamped P a p if 

2cos(2A) > 1 + k. (31) 

Note that if tan(#) = 1, then Eqs. ()29|) and (}3~T]) will have the same form except for the sign 
of the left-hand side. 

In Fig. Q the qualitative effects of neutrino wave packet decoherence and neutrino decay on 
the neutrino survival probability are shown. From this figure, we clearly see how the wave 
packet decoherence simply corresponds to a damping of the oscillating term and the decay 
of all mass eigenstates corresponds to an overall damping of the undamped neutrino survival 
probability. For the case of only one decaying mass eigenstate, the probability converges 
towards the square of the content of the stable mass eigenstate in the initial neutrino flavor 
eigenstate. 

3.2 Three-flavor electron- muon neutrino transitions 

For a fixed neutrino oscillation channel, the damped neutrino oscillation probability Eq. (jSJ) 
can be written more explicitly in terms of the mixing parameters and the mass squared 
differences. Below, we will use the standard notation for the leptonic mixing angles, i.e., 



13 




2 4 6 8 



Oscillation phase (A) 



Figure 1: The qualitative effect of different damping signatures on the two- flavor neutrino survival 
probability as a function of the oscillation phase A. The mixing used in this plot is maximal (9 — n/2) and 
the damping parameters have been highly exaggerated. The scenario "Oscillation + decay I" corresponds 
to decay of both mass eigenstates with equal rates, whereas "Oscillation + decay II" corresponds to the 
second mass eigenstate decaying while the first mass eigenstate is stable. 



Sij = sin(#jj) and Cy = cos(%). Then, for example, the v e survival probability P ee is given 
by 



P - r 4 

1 ee — C 13 



1 



D u c 4 12 + D 22 s\ 2 + ]-D 21 sin 2 (2# 12 ) cos(2A 21 ) 



+ ^ sin 2 (2# 13 ) [D 31 c\ 2 cos(2A 31 ) + D 32 s\ 2 cos(2A 32 )] + D 33 s$ 3 , 



(32) 



which is dependent on all neutrino oscillation parameters except for 6 23 and Sep, while the 
probability P eil of oscillations into is given by 



P 



e/i 



- sin 2 (2^ 12 )c 2 3 [( J D 11 + D 22 ) - 2D 21 cos(2A 21 )] 

+ i sin(2# 12 ) sm(29 23 ){c s [D llC 2 12 - D 22 s\ 2 - D 21 cos(2fl 12 ) cos(2A 21 )] 

-D 21 s s sin(2A 2 i) + D 32 cos(2A 32 - 5 C p) - Asi cos(2A 3 i - 5 C p)} S13 
+s 2 23 [D u c 4 12 + D 22 s\ 2 + D 33 - 2D 3l s\ 2 cos(2A 31 ) - 2,D 32 c 2 2 cos(2A 32 )] s 2 3 

+ ^ sin 2 (2^ 12 )[2D 21 cos(2A 21 ) - cj 3 (D u + D 22 )\ s\ 3 + 0(s 3 13 ), 



(33) 
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where sg = sin(<5cp) and c§ = cos(5cp). Furthermore, the survival probability can be 
computed to be of the form 



- sin 2 (2# 23 ) [D 32 c\ 2 cos(2A 32 ) + D 31 s 2 12 cos(2A 31 )] 



+ c 23 



D ix s{ 2 



+ D 22 c{ 2 + -£> 21 sin 2 (2# 12 ) cos(2A 21 ; 



^33S 23 



+c s sin(2^i 2 ) sin(2# 23 ) {c 2 23 [D n s 2 12 - D 22 c 2 12 + D 21 cos(26 12 ) cos(2A 21 ; 
+s 2 23 [D 31 cos(2A 31 ) - D 32 cos(2A 32 )] } s 13 + 0{s\ 3 ). (34) 



Note that the probabilities P efl and P^ are series expansions in s 13 , whereas the probability 
P ee is valid to all orders in s 13 . The reason to use these expressions rather than the exact 
expressions is that, unless some further assumptions are made, the formulas for P e/X and P w 
are quite cumbersome. 



The probability P Me can be obtained by making the transformation 5c p 
probability P t 



fie 

' efJi , i.e., P Me = P e fj,(Scp — > Sep)- Furthermore, in vacuum, the anti-neutrino 
oscillation probabilities can be obtain from the neutrino oscillation probabilities through 
the same transformation as above. Note that this is not true for neutrinos propagating in 
matter. 



-Sep in the 



3.3 Probabilities for decoherence-like effects in experiments 



For a decoherence-like damping effect, Da = 1 for all % and the relations 

£ P al 3 = l and J2 P ^ = 1 ( 35 ) 

ct=e,fi,T f3=e,/j,,T 

are still valid despite the presence of damping factors (i.e., no neutrinos are lost due to 
effects such as invisible decay, absorption, etc.). Note that, in the case of a decoherence-like 
damping effect, all neutrino oscillation probabilities can be constructed from P ee , P e[1 , and 
Ppp due to the conservation of total probability given in Eq. (|33j) . 

It is interesting to observe what effect a decoherence-like damping could have on the neutrino 
oscillation probabilities for different experiments. Therefore, we will now study different 
kinds of neutrino oscillation experiments and make different approximations depending on 
the type of experiment to investigate what the main damping effects are. 



Short-baseline reactor experiments 

Short-baseline experiments, such as CHOOZ [66,67] and Double-CHOOZ [68], are operated 
at the atmospheric oscillation maximum A 3 i ~ A 32 = 0(1) in order to be sensitive to 
sin 2 (2#i 3 ). The most interesting quantity is the v e survival probability P ss . For these 
experiments, it turns out (see Sec. 0]) that it is important to keep all damping factors. As 
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a result, the v e survival probability is given by 

Pee = ^sin 2 (2M[l-^2icos(2A 21 )]| 

+ ^sm 2 (2^3)[ J D3iC? 2 cos(2A 31 ) + J D 32 ^ 2 cos(2A 32 )] + 4- (36) 

The most apparent feature of this equation is the term within the curly brackets, which has 
the form of the survival probability for a two-flavor neutrino damping scenario with 9 = 9\ 2 
and A = A 21 . Therefore, even in the limit #13 — > [close to the sin 2 (26 l 1 3) sensitivity limit], 
the damping factor D 2 i might be constrained by the contribution of the solar oscillation at 
low energies. Furthermore, in the limit A 2 i — > (or large #13), D 2 % is close to unity [cf., 
Eq. (0)] and D 3 % ~ D 32 (this could be expected if A21/A31 — > 0), then this expression will 
exactly mimic the two-flavor neutrino damping scenario with 9 = # 13 and A = A 31 = A 32 . 
Thus, depending on which small number (the ratio of the mass squared differences or si 3 ) 
is the largest, two different two-flavor neutrino scenarios are obtained as expected from 
the non-damped case. If 9\ 3 is relatively large (compared to the ratio of the mass squared 
differences), then the latter two-flavor case will apply. It is then interesting to note that 
the damping factor D 3 i, the neutrino source energy spectrum, and the cross-sections all 
have some energy dependence, which means that they can "emphasize" certain regions in 
the energy spectrum which are most sensitive to damping effects. If we assume that the 
total impact is strongest close to the oscillation maximum, then the damping effect will be 
misinterpreted as a smaller value of sin 2 (26*i 3 ) [cf., Eq. ([21)1. which will in both cases be 
closer to unity]. Therefore, as we will demonstrate, any such damping can fake a value of 
sin 2 (26 l 13 ) which is smaller than the one that is provided by Nature. 

Note that, for the case of wave packet decoherence, D 21 , D 32 , and D 31 are not independent 
[cf., Eq. (J7|)], which means that any of the terms in Eq. (|36|) could lead to information on 
the parameter oe- 



Long-baseline reactor experiments 

For long-baseline reactor experiments operated at the solar oscillation maximum A 21 = 
(9(1), such as the KamLAND experiment [5,6], the damping factors D 3 i and D 32 of a 
decoherence-like scenario with £ > are small, since the large mass squared difference 
makes the argument of the exponential functions in Eq. (jlj) large and negative. In addition, 
these two damping factors are attached to neutrino oscillations associated with the large 
phases A 31 and A 32 [see Eqs. ()32jl -(jM j) ]. which effectively average out. As a result of these 
two effects, the oscillating terms involving the third mass eigenstate can be safely set to 
zero. After some simplifications, the v e survival probability P Se - is found to be 

Pee = 4 jl - i sin 2 (2# 12 )[l - D 21 cos(2A 21 )] j + s 4 13 . (37) 

This expression is clearly of the familiar form P g g = cf 3 P~ + sf 3 , where P|| is the damped 
two-flavor v e survival probability with 9 = 9\ 2 and A = A 21 , which is also obtained in the 
non-damped case when averaging over the fast oscillations [cf. Eq. (|36|)]. For the case of wave 
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packet decoherence, we know from Eq. (JJJ) that the parameter oe could be constrained by 
either of these two equations. Since this parameter is experiment dependent, one could argue 
that one should obtain some limits from the KamLAND experiment, because the reactor 
experiments are very similar in source and detector (see, e.g., Ref. [69]). However, it should 
be noted that KamLAND has a rather weak precision on the corresponding 9 12 measurement 
because of normalization uncertainties. Since a decoherence contribution would appear at 
low energies, the data set in Ref. [6] does not seem to be very restrictive for the parameter 
a E . 

Beam experiments 

For beam experiments, such as superbeams, beta-beams or neutrino factories, one may 
assume A21 — as a first approximation if one wants to be sensitive to sin 2 (2^13) , since, 
at the energies and baseline lengths involved, the low-frequency neutrino oscillations do 
not have enough time to evolve. In the case of £ > 0, this also implies that _D 12 = 1 and 
D = D 32 = D31 to a good approximation. From these assumptions, it follows that 

P e » = 24[l- J Dcos(2A)] S 2 3 + 0(4), (38) 

P w = l-isin 2 (2^ 3 )[l-^cos(2A)]+0( S 2 3 ), (39) 

where A = A 32 = A31. Note that the probability P etl is correct up to 0(sf 3 ) [as compared 
with Eq. (J3l?j) . which is only valid up to 0(s 2 3 )}, this is one of the cases where the assumptions 
made simplifies the s 2 3 term in this probability. Both of the above equations show obvious 
similarities with the cases of damped two-flavor neutrino oscillations. For P efl we have an 
approximate two-flavor neutrino scenario with s 2 c 2 = s 23 sf 3 and P w is a pure two-flavor 
neutrino formula with 9 = # 23 up to the corrections of order s 2 3 . Since the disappearance 
channel P m at a beam experiment is supposed to have extremely good statistics, D will be 
strongly constrained by this channel. Note that the damping in P w qualitatively behaves 
as the one in Eq. (}2~T]) . i.e., the damped probability might be larger or smaller than the 
undamped probability depending on the position relative to the oscillation maximum E max . 

3.4 Probabilities for decay-like effects in experiments 

If £ = and an 7^ 0, then Da 7^ 1 and Eq. (|3~K|) will not hold. We define any effect of this 
kind to be "probability violating". As mentioned in the two-flavor neutrino discussion, a 
very interesting special case of the probability violating effects is the case of a decay-like 
effect. The neutrino oscillation probabilities for decay-like effects corresponding to the ones 
given for decoherence-like effects are listed below. 

Short-baseline reactor experiments 

For the short-baseline reactor experiments, we obtain the v e survival probability as 

P s -e = c\ 3 {{A lC \ 2 + A 2 s\ 2 ) 2 - A l A 2 sin 2 {29 l2 ) sin 2 {A 2 i)} 

+A 3 s 2 13 {A 3 s 2 13 + 2c 2 13 [A lC 2 12 cos(2A 31 ) + A 2 sj 2 cos(2A 32 )]}. (40) 
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Again, as in the case of decoherence-like damping, the expression within the curly brackets 
is of a two-flavor form with 6 = 9 12 and A = A 12 . In the limit when sin 2 (26 , 13 ) is large and 
we ignore the solar oscillations, we obtain the two-flavor neutrino scenario 

p.. = A 2 {(c 2 3 + ks 2 13 ) 2 - Ksin 2 (2# 13 ) sin 2 (2A)} (41) 

only if we assume that A\ = A 2 = A, where A = A31 = A32 and k = A3 /A. 

Long-baseline reactor experiments 

Assuming that the fast neutrino oscillations average out, the v e survival probability is given 
by 

P-e-e = 4^1 + AMs, (42) 

where P~ is the two-flavor decay-like v e survival probability with 9 = # 12 and A = A 2 i 
[cf., Eq. (J37|) ]. In this expression, the sf 3 term is also damped, which does not apply in a 
decoherence-like scenario. 

Beam experiments 

When the assumptions A21 — and A = A\ = A 2 (which could be expected in a decay 
scenario where m\ = m 2 ) are made, the neutrino oscillation probabilities that are relevant 
for beam experiments become 

P eM = A 2 s 2 23 [l + k 2 ~2Kcos(2A)]s 2 13 + 0(s :i 13 ), (43) 
= A 2 [(c 2 23 + Ks 2 23 ) 2 -Ksm 2 (2e 23 )sm 2 (A)]+0(si 3 ), (44) 

where k = A 3 j A and A = A 32 = A 31 . These probabilities mimic decay-like two- flavor 
probabilities just as the corresponding decoherence-like effects mimic decoherence-like two- 
flavor probabilities to leading order in S13. 

4 Application I: Faking a small sin 2 (20i 3 ) at reactor experiments 
by decoherence-like effects 

In this section, we demonstrate the possible effects of damping at a simple example using a 
full numerical simulation. Let us only consider the case of intrinsic wave packet decoherence, 
which is very interesting from the point of view that it is a "standard" effect in any realistic 
neutrino oscillation treatment. However, similar effects could occur from related signatures, 
such as quantum decoherence. As experiments, one could, in principle, consider all classes of 
experiments in order to investigate decoherence signals. New reactor experiments with near 
and far detectors [70,71] are candidates for "clean" measurements of sin 2 (26 l 1 3), i.e., they 
are specifically designed to search for a sin 2 (26*13) signal. As we have discussed in Sec. 13.31 
an interesting decoherence-like effect at such an experiment would be a derived value of 
sin 2 (2^i3) which is smaller than the value provided by Nature. In this case, the CHOOZ 
bound might actually be too strong and the interpretation of new reactor experiments might 
be wrong. 
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Figure 2: The neutrino oscillation probability P g g (left) and event rates (right) for the experiment 
Reactor-I from Ref. [71] in the analysis range. For the simulated parameter values, we use Am 2 ^ = 
2.5- lCT 3 eV 2 , Am|i = 8.2- lCT 5 eV 2 , sin 2 29 12 = 0.83, sin 2 29 23 = 1, S C p = [22,72-74] and the values for 
sin 2 (26 l i3) and <je as given in the left plot. 

If we assume that there is an intrinsic loss of coherence, then the reactor v e survival proba- 
bility Pgg will be given by Eq. (|36|). In order to illustrate the decoherence effect, we show in 
Fig. El -Pee and the corresponding event rates for the experiment Reactor-I from Ref. [71] 
(full analysis range shown). The different curves correspond to the non-oscillatory case 
as well as different combinations of sin 2 (2# 13 ) and <je- As one can observe, the two cases 
sin 2 (26 , 13 ) large and decoherence [sin 2 (26*i 3 ) = 0.05 and oe = 2MeV] and sin 2 (20i 3 ) small 
and no decoherence [sin 2 (26 ) i 3 ) = 0.03 and <je = 0] correspond, especially in the event rate 
plot, very well to each other [as compared to the other two cases of no oscillations and 
large sin 2 (2^13) only]. This means that the decoherence effect can mimic a smaller value 
of sin 2 (26 , i 3 ) than what is provided by Nature. Note that in the probability plot, there is 
a significant contribution from loss of coherence in the solar terms for low energies. As we 
will see later, this contribution can limit the decoherence effects even for no sin 2 (26 l i3) sig- 
nal. In addition, the damped neutrino oscillation probability is larger than the undamped 
one in the range discussed after Eq. (|21|). where the oscillation maximum is here at about 
£ max ~3.4MeV. 

In order to illustrate the effect for a complete analysis, we show in Fig. |3] the simultane- 
ous sensitivity to sin 2 (20 13 ) and o E for Reactor-I (£ = 400tGWyr) and Reactor-II 
(C — 8 000 1 GWyr) from Ref. [71] (1 d.o.f.) using an extended version of the GLoBES soft- 
ware [77]. In this figure, ge is assumed to be a free (fit) parameter that has to be measured 
by the experiment. Therefore, without additional knowledge, the sin 2 (2#i 3 ) sensitivity limit 
is obtained as a projection of the curves onto the sin 2 (2#i 3 )-axis. Since the sin 2 (2#i3) sensi- 
tivity limit for no decoherence effects is the one for ge — 0, the arrows indicate the shift of 
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Figure 3: Simultaneous sensitivity to sin 2 (2#i3) and o\e for the experiments Reactor-I (left) and 
Reactor-II (right) from Ref. [71] (curves shown for 1 d.o.f). For the simulated parameter values, we use 
sin 2 (26'i3) = 0, oe = 0, and the other values as in Fig. |2 For the thick solid curves, the unshown fit 
parameter values are marginalized over, where post-KamLAND external precisions of 10 % on each Am^i 
and 9\2 [75,76] are imposed along with an external error of 10 % for Am^ obtained from the superbeams 
is assumed. For the thin dashed curves, the unshown fit parameter values are fixed (no correlations). For 
the numerical analysis, an extended version of the GLoBES software [77] is used. The arrows indicate the 
shift of the sin 2 (26*i 3 ) sensitivity limit if one assumes oe as a free parameter. 

this limit by the unknown oe- This means, for example, that the sensitivity limit becomes 
about 50 % to 100 % worse than that for the actual oe = 0, since the decoherence mimics a 
smaller value of sin 2 (2^13) than what is provided by Nature. Similar results to the left plot 
are obtained for the proposed Double-CHOOZ experiment [68]. Note that the correlation 
between sin 2 (2# 13 ) and oe affects the sin 2 (2# 13 ) sensitivity (projection onto the horizontal 
axis), but not the oe sensitivity (projection onto the vertical axis). The latter is correlated 
with the other neutrino oscillation parameters (especially the solar parameters), as one can 
read off from the difference between the solid and dashed curves. For the <je sensitivity, one 
obtains a E < 10 MeV (Reactor-I) and a E < 5 MeV (Reactor-II) at the 3a confidence 
level. As one can observe from the left plot of Fig. |2J there is some contribution of the 
solar oscillation averaging to the decoherence effect at low energies. In fact, this is the 
reason why one can constrain cte even for sin 2 (2#i3) = 0, since in the decoherence effect, the 
atmospheric oscillations are suppressed by the oscillation amplitude sin 2 (26 \z) . Obviously, 
this solar decoherence effect determines the upper bound for cr^, which means that the <te 
sensitivity is limited by the knowledge on the solar oscillation parameters [cf., Eq. (I36|) ]. 

As we have discussed in Sec. |2J <Je might be an experiment dependent parameter related 
to the production and detection processes. Instead of deriving bounds for this parameter 
from reactor experiments, one can estimate from Fig. El that one has to constrain <je better 
than to about oe < 0.5 MeV in order not to have a significant deterioration of the sin 2 (2#i3) 
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Figure 4: Contributions of the first three different damping signatures from Tabled to the disappear- 
ance probability as function of the neutrino energy. Here L — 3 000 km and the neutrino oscillation 
parameters as in Fig. [21 with sin 2 (26*13) = are used. The parameters for the non-standard effects are given 
in the plots, where zero corresponds to the thick curves (oscillations only) and larger values correspond to 
curves further off the zero curve. The energy range corresponds to the analysis range of the 50 GeV neutrino 
factory NuFact-II from Ref. [78]. 

sensitivity limit. In addition, in order to exclude an experiment dependent effect, it is highly 
recommendable to measure the same quantity with different techniques such as sin 2 (2(9i 3 ) 
with both reactor experiments and superbeams. 



5 Application II: Testing and disentangling damping signatures 
at neutrino factories 

If we want to constrain the model parameters in Table ^ and to test the different models 
against each other, then we will need to choose a high-precision instrument to test these 
tiny effects. Therefore, we investigate the potential of a neutrino factory. In particular, the 
muon neutrino disappearance channel — > at a neutrino factory has very good statistics 
and the impact of neutrino oscillation parameter correlations other than with Am^ and 
#23 is very small. Thus, we will mainly focus on this disappearance channel, but include 
the appearance information in the full analysis and demonstrate how the value of sin 2 (2#i3) 
would influence the effects. Since our exponential damping model is not directly comparable 
to other approaches in the literature, we put a major emphasis on the identification problem 
of a non-standard contribution: If we actually observe something unexpected, how well can 
we determine what sort of effect this actually is? In the simplest case, this means that 
we test a signature against the standard (no damping) scenario giving us limits for the 
parameters. Since it is almost impossible to include the correlations among all parameters, 
we choose to use = a independent of i and j in this section in order to drastically reduce 
the number of parameters. This means that we now have to deal with eight correlated 
parameters (six neutrino oscillation parameters, the matter density, and the parameter a). 
We have motivated this choice at the end of Sec. 12. II and, for individual cases, in Sec. 12.31 

Before we come to the results of a complete simulation, let us illustrate the spectral behavior 
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Simulated damping signature 


Fit signature 


Decoherence Decay Oscillations Absorption Q. decoh. I Q. decoh. II 

GcV ~ in-5 GcV ~ in- 7 nV 4 ~ lo- 8 ~ lo-i" ^ io^ 4 ~ 

km GeVkm GcV^km ~aV T 


No damping 


1.7 (2.8) 4.3 (7.2) 5.1 (8.3) 1.9 (3.1) 4.1 (6.8) 2.0 (3.6) 


Decoherence 
Decay 
Oscillations 
Absorption 
Q. decoli. I 
Q. decoh. II 


4.3 (7.2) 5.1 (8.3) 1.9 (3.1) 4.1 (6.8) 2.0 (3.6) 
1.7 (2.8) - 6.3 (10) 3.4 (5.7) 6.0 (10) 2.6 (5.1) 
1.7 (2.8) 5.8 (9.8) - 1.9 (3.2) 4.1 (6.9) 13 (17) 
1.7 (2.8) 7.8 (13) 5.2 (8.5) - 24 (40) 2.1 (3.8) 
1.7 (2.8) 6.3 (11) 5.1 (8.3) 11 (19) - 2.1 (3.7) 
1.7 (2.8) 4.3 (7.2) 5.1 (8.3) 1.9 (3.1) 4.1 (6.8) 


All models 


1.7 (2.8) 7.8 (13) 6.3 (10) 11 (19) 24 (40) 13 (17) 



Table 3: Parameter sensitivity limits for which the simulated models (in columns) from Table^could be 
distinguished from the fit models (in rows) at the 3cr (5c) confidence level (for the experiment simulation 
NuFact-II from Ref. [78]). For example, decoherence could be established against all models (including 
standard oscillations) for the simulated > 1.7 GeV. For the simulated neutrino oscillation parameter 
values, we use the same values as in Fig. Hand sin 2 (26»i 3 ) = as given in the column captions. The 
fit parameter values (including the model parameter a) are marginalized over. The row "no damping" 
corresponds to the standard neutrino oscillation scenario, i.e., it corresponds to the upper bounds for the 
parameters assuming that there is only one non-standard effect. The row "All models" corresponds to the 
most conservative case, i.e., it is an estimate for how well one can establish the model against all of the 
other shown models. 

(energy dependence) of the neutrino oscillation probability P w in Fig. 0] for some charac- 
teristic examples. Earlier in Sec. El we have already discussed that there are two general 
interesting cases: Either only the oscillatory terms are damped or all terms are damped. 
In Fig. HJ we can clearly identify this difference between the decoherence-like and the other 
two damping models (decay and oscillations). In all the shown cases (for which 7 > 0), 
the relative importance of the damping increases as the energy decreases. However, since 
also the neutrino oscillation probability drops with lower energies, the absolute size of the 
effect is determined by the ratio of signature versus probability effect for low energies. In 
addition, cross-section and flux will disfavor low energies, which means that the low-energy 
effects become even harder to identify. This makes the wave packet decoherence scenario 
most difficult to test, since the E~ 4 dependence in the exponent strongly favors low ener- 
gies. However, it might be most easily distinguished from the decay and oscillation damping 
scenarios because of its unique signature. As we have discussed after Eq. (|2Tj) [which also 
holds for the similar Eq. (j39j) ] . it is a characteristic feature of decoherence-like signatures 
that they cross the undamped curve at 2E max /3 and 2E max , which here evaluate to 4 GeV 
(outside of the analysis range) and 12 GeV. In Fig. 0] (left panel), this effect is hardly ob- 
servable because of the E~ A energy dependence, but the quantum decoherence motivated 
case "Quantum decoherence II" from Table clearly shows this behavior because of an E~ 2 
energy dependence. As far as the other two signatures are concerned, the decay damping 
has a linear energy dependence in the exponent as opposed to the quadratic one for the 
oscillation damping scenario. Therefore, one has the strongest high-energy effect for the 
decay damping scenario. 

In order to test the different models against each other, we use a modified version of the 
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GLoBES software [77] and the neutrino factory setup NuFact-II from Ref. [78]. This neu- 
trino factory uses a 50 kt magnetized iron detector, 4 year of running time in each polarity, 
and 4MW target power (corresponding to 5.3 • 10 20 useful muon decays per year). For a 
fixed set of simulated parameter values including the simulated damping parameter a, we 
marginalize over the fit neutrino oscillation parameters including the fit damping parameter. 
Due to the complexity of the parameter space, we assume that the sgn(Am| 1 )-degeneracy 
has been resolved by this time. We define the sensitivity limit to a as the threshold above 
which the simulated damping model could be distinguished from the fit damping model. 
Thus, if the damping mechanism is really there, then the damping parameter a has to be 
above this threshold in order to establish the model against the fit model with the given 
experiment. In particular, we include the fit damping model "no damping", which corre- 
sponds to the standard neutrino oscillation case. For the simulation, we impose external 
precisions of 10 % on each 12 and Am^ [75,76]. In addition, we assume a constant matter 
density profile with 5 % uncertainty, which takes into account matter density uncertainties 
as well as matter density profile effects [79-81]. However, we assume that the neutrino 
factory itself measures Am^ and #23 with its disappearance channel, i.e., we do not impose 
an external precision on these parameters. 

The resulting sensitivity limits of this analysis are shown in Table where the columns 
correspond to the simulated models and the rows correspond to the fit models. These 
results are computed for sin 2 (26 , 13 ) = 0. It turns out that for a simulated value of sin 2 (26 ) i 3 ) 
close to the CHOOZ bound sin 2 (26> 13 ) ~ 0.1, the limits on a would improve up to about 
30 % [depending on model and value of sin 2 (26'i3)] because of the additional contribution 
from the appearance signal. 5 Let us first of all discuss the resulting sensitivities against the 
standard neutrino oscillation scenario for some simple cases. For decoherence, the obtained 
numbers indeed correspond very well to the energy resolution of the detector, which is about 
15 % of the neutrino energy, i.e., 1.5 GeV for a neutrino energy of E = 10 GeV, where the 
major effect takes place (c/., Fig. left plot). Since the neutrino oscillation probability 
changes sufficiently fast in this region, the measurement is limited by the energy resolution 
of the detector. In the wave packet approach, the bound against the "no damping" model 
a E ^ 1.67 GeV translates into a x > 6- 10~ 17 m. This rather small number (sub-nucleon size) 
means that the bound is not very useful for wave packet decoherence, since it is virtually 
impossible to create such sharply peaked wave packets. However, there might be other 
energy averaging effects that can be constrained. For decay, we obtain a limit, against the 
standard model, which is comparable to the current neutrino lifetime limit for m^. Note 
that we have included all correlations with the neutrino oscillation parameters in this limit. 
However, the limit would be a factor of two weaker if we considered only decay of 777.3 instead 
of all mass eigenstates. Since there are quite strong bounds on the mi and m 2 lifetimes 
from supernova and solar neutrino observations, this factor of two difference should be a 
very good approximation for the actual limit. For the oscillation signature, the obtained 
limits are of the order of magnitude 5 • 10~ 7 eV 2 , which corresponds to (Am 2 3 ) 2 times the 
active-sterile mixing in our estimate for a possible mass scheme (c/., Sec. 12. 3|) . Considering 
the Am 2 ^ dependence, this is in fact not a very strong bound. However, note that we 

5 We do not show these results, since the exact interpretation of the appearance signal is model dependent. 
In addition, matter effects are strong in this case and they depend on the treatment of those in the context 
of the damping model. 
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Correlations among different effects 
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Figure 5: The impact of different correlations on the statistics (and systematics) sensi- 
tivity limit of the model dependent parameter a (3c), where the horizontal axis represents 
multiples of the statistics (and systematics) sensitivity limit. The group captions refer to 
the simulated models and the bar labels to the fit models, where only the fit models are 
shown which affect the sensitivity limit more than by 5 %. The dark bars represent the 
correlations with the neutrino oscillation parameters (fit parameter a = fixed) and the 
light bars indicate the additional change if the model specific fit parameter a is marginalized 
over. The lowest light bar extends to 37. 

have taken into account the full parameter correlation, i.e., this effect could not come from 
sin 2 (2^13) or any other standard parameter. 

In order to discuss the general identification problem among different damping signatures, 
some information can be obtained from Table El In addition, in Fig. El we show the impact 
of the correlations with the standard neutrino oscillation parameters (dark bars) as well as 
the additional correlation with the fit model parameter a (light bars) on the a sensitivity 
limit for the simulated models from Table ^ The horizontal axis shows the ratio of the 
a sensitivity limit including correlations to the one from statistics and systematics only 
(which corresponds to 1), where we only include fit models with relevant model parameter 
contributions. Two models are highly correlated if a possible signature in one model can 
be compensated by a change of parameter(s) in the other. Since we include the standard 
neutrino oscillation scenario in all models, a small change in the fit neutrino oscillation 
parameters might also compensate a damping signature within the measurement precision 
of the experiment. Therefore, we include for all signatures the standard neutrino oscillation 
parameter correlation as dark bars, i.e., the dark bars represent the fit against the standard 
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neutrino oscillation scenario (for fixed fit parameter a), and the light bars are a measure 
for the additional problem to distinguish a non-standard signature from the ones of other 
possible non-standard models. The interpretation of these bars is as follows: The dark bars 
reflect the limit (right edges) for a (as multiple of the statistics limit) beyond which the 
non-standard signature could be distinguished from the standard neutrino oscillation case 
at the 3a confidence level. However, if a should be within one of the light bar ranges, then 
it could not be uniquely identified, since it could also well be the non-standard signature 
corresponding to this bar. From Fig. we make a number of interesting observations: 

• Signatures which have negative 7's (Absorption and Quantum decoherence I) are 
almost not affected by correlations with the neutrino oscillation parameters, i.e., they 
cannot be explained by different neutrino oscillation parameter values. In these cases, 
the spectrum is more suppressed for large values of E than for small values, which 
means that the signature behaves unlike an oscillation signature corresponding to 
7 = 2. However, it is difficult to identify which of these models is realized. 

• Signatures with 7 = 2 (Oscillations into v s and Quantum decoherence II) are highly 
affected by correlations with the standard neutrino oscillation parameters, since the 
signatures have an energy dependence similar to the oscillation signature. Similar 
signatures, such as decay, can enhance this correlation. 

• Unique signatures (Wave packet decoherence and Neutrino decay) can easily be dis- 
tinguished from all the other models. Although there could be some correlations with 
similar signatures for neutrino decay, the absolute impact on the a sensitivity limit is 
comparatively small (up to a factor of three). 

6 Summary and conclusions 

We have introduced exponential damping factors in the neutrino oscillation probabilities, 
which lead to distinctive signatures, i.e., energy dependent damping effects in the energy 
spectrum. These damping factors are one approach to test non-oscillation effects on the 
neutrino oscillation probability level. They can be motivated by many different models 
such as intrinsic wave packet decoherence, neutrino decay, oscillations into sterile neutrinos, 
neutrino absorption, quantum decoherence, etc.. They describe the second order contribu- 
tions of small possible "non-standard" corrections to the three-flavor neutrino oscillation 
framework (in vacuum as well as in matter) on a rather abstract level. As opposed to 
tests of probability conservation, the damping factors can, in addition, describe a damping 
of the oscillating terms (which preserves the total probability) as well as they imply, by 
their energy dependence, some information on the type of effect. We have demonstrated 
how damping factors can modify the neutrino oscillation probabilities relevant for future 
high-precision short- and long-baseline experiments, since these experiments might be most 
sensitive to very small spectral effects. 

As one application, we have shown that decoherence- like damping signatures can severely 
modify the interpretation of experiments, where we have chosen wave packet decoherence 
damping at new short-baseline reactor experiments as an example. In this case, two com- 
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peting small effects, namely the effect of a non-zero value of sin 2 (26*13) and a damping con- 
tribution, might be mixed up. In particular, the damping could fake a value of sin 2 (26> 13 ) 
which is much smaller than the value provided by Nature. Such a sin 2 (26> 13 ) suppression 
effect can either be intrinsic (such as quantum decoherence) , experiment dependent (such 
as some averaging effect not taken into account), or both (such as wave packet decoherence 
related to the production and detection processes). Intrinsic effects will be observable by 
all types of experiments, which means that there are very stringent limits available from 
existing data as well as future experiments will test the consistency of the picture. On the 
other hand, experiment dependent effects can only be checked by complementary techniques 
measuring the same quantity. One such complementary pair has, in the past, been the solar 
and long-baseline reactor experiments. In the future, it will therefore be very important to 
measure sin 2 (26> 13 ) by reactor experiments and superbeams as complementary techniques, 
since one of them alone could fail for such experiment dependent effects. Eventually, the 
LSND experiment could be a strong hint for such an experiment dependent effect if it is 
rejected by the MiniBooNE experiment. 

One of the most interesting features of damping signatures are their characteristic spectral 
(energy) dependencies, which can act as a "fingerprint" for many sources of non-oscillation 
effects. For example, specific signatures could point to new interesting physics beyond the 
standard model. We have therefore discussed how large the effects from different damping 
signatures have to be in order to be identified and how well these damping signatures 
could be distinguished for the example of neutrino factories. In some cases, such damping 
signatures can be compensated by a shift of the neutrino oscillation parameters, which means 
that given such a damping effect, it is quite likely to obtain an erroneous determination of 
these parameters. However, if the damping effects are strong enough, then an establishment 
of non-oscillation effects will be possible. Once such a damping effect is established, it will 
be very interesting to know from which non-standard mechanism it actually arises. Given 
this question of the identification problem, we have found that signatures with a damping 
similar to exp(— aL 13 /E 1 ), 7 = 1, 2, . . . are strongly correlated (peaking at 7 = 2) with the 
standard neutrino oscillation parameters, i.e., it is difficult to distinguish them from small 
adjustments in the neutrino oscillation parameters. However, damping signatures similar to 
exp(— aL^E 2 ) can be very easily disentangled from the neutrino oscillation parameters, but 
it is difficult to distinguish them from each other. It is also extremely difficult to establish 
a damping of the oscillations against a damping of the probabilities with the same spectral 
index 7 because of the correlations with the neutrino oscillation parameters. 

Finally, we conclude that spectral tests of damping signatures in neutrino oscillation prob- 
abilities are an important test of the consistency of the three-flavor neutrino oscillation 
picture. If any deviation from this picture is found, then the most important question will 
be what sort of effect we are dealing with. Exactly this information could be provided by 
the spectral dependence of the damping signature, which means that this approach could 
be an important test of physics beyond the standard model. 
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